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Abstract 

We investigate the phase diagram of strongly interacting matter as a function of temperature 
and baryonic density /chemical potential, within Nambu-Jona-Lasinio type models. We perform a 
systematic study concerning the existence, location, and properties of a critical end point /tricritical 
point, both in SU(2) and SU(3) versions of the model. We verify that, for = md = and up 
to a critical strange quark mass, there is a tricritical point, which becomes a critical end point in 
a world with realistic values of the current quark masses. The properties of physical observables, 
such as the baryon number susceptibility and the specific heat, are analyzed in the vicinity of 
the critical end point, with special focus on their critical exponents. The behavior of mesons in 
the T — ^b{pb) plane is analyzed in connection with possible signatures of partial and effective 
restoration of chiral symmetry. 
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I. INTRODUCTION 



Recently there has been encouraging progress on nonperturbative studies of the QCD 
thermodynamics which have stimulated a great deal of theoretical activity. Phenomenolog- 
ical and microscopic models have been developed along parallel and complementary lines 
allowing to predict a rich phase structure at finite temperature, T, and chemical potential, 



The quark gluon plasm (QGP) is a longstanding theoretical issue since the 
discovery of the asymptotic freedom of QCD ^, 5]. Besides the intrinsic theoretical inter- 
est of this subject, such studies are important because they are directly applicable to the 
regime under current experimental investigation at the Brookhaven National Laboratory 
(BNL) Relativistic Heavy Ion Collider (RHIC). In fact, extensive experimental work has 
been done with heavy-ion collisions at CERN and Brookhaven to explore the T ~ fiB phase 
diagram and look for signatures of the QGP. 

Theoretical studies have been accumulating a lot of evidence that there exists a critical 
end point (CEP) in the phase diagram of strongly interacting matter. Since Fodor and Katz, 
who presented a phase diagram with the CEP within lattice calculations 6|, remarkable 
progress in this area has been made. It is an open question, whether a critical end point 
exists on the T — fis plane and, particularly, how to predict its location. When only thermal 



and lattice simulations 



9| indicate that the 



effects are concerned, universal arguments P, Q] 
order of the phase transition depends on the masses and flavors of quarks. 

Considering also nonvanishing chemical potentials, a variety of models (see e.g. 
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III) 



predict a second order phase transition point in the phase diagram. This suggests that the 
phase diagram exhibits a CEP. At this point the phase transition is of second order and long 
wavelength fluctuations appear, leading to characteristic experimental consequences that 
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13j. So, the 



can be detected by enhanced critical fluctuations in heavy-ion collisions 
location of the CEP has become an important topic in effective model studies and lattice 
calculations. In fact, the phase diagram and QCD thermodynamics in general are becoming 
more transparent due to the combination of research in several areas: perturbative QCD, 
effective models, and lattice calculations. 

The possible existence of such a point has recently been emphasized and its universal 
critical properties have been discussed by several authors in the context of QCD inspired 
models l3. 
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14 



15| . This point of the phase diagram is the special focus of the present 
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article. 

In a previous work [15|], we studied the phase diagram focusing our attention on the 
CEP and the physics near it, through the behavior of the baryon number susceptibihty 
and the specific heat; the study was performed in the framework of the SU(3) Nambu- 
Jona-Lasinio (NJL) modeL Here, besides extending the investigation to other observables, 
we make a comparative study of the phase diagram in the SU(2) and SU(3) NJL models. 
Since more information can be taken within the simpler version of the NJL model, this 
systematic study is expected to provide a better understanding of the interesting physics 
around the CEP/TCP (tricritical point). Our main goal is to locate the critical end point 
and confront the results with universality arguments. Based on the fact that the CEP 
is a genuine thermodynamic singularity, being considered a second order critical point, the 
order parameter and related observables, like susceptibilities, can provide relevant signatures 
for phase transitions. We notice that susceptibilities in general are related to fluctuations 
through the fluctuation dissipation theorem, allowing to observe signals of phase transitions 
in heavy-ion reactions IgI . The specific heat C, which is related to the event-by-event 



temperature fluctuation [17l], and mean transverse momentum fluctuations [18|] in heavy- ion 
reactions, is also a quantity of interest in our calculation. These fluctuations should show a 
divergent behavior near the CEP. After equilibration, the dense matter created in relativistic 
heavy-ion collision will expand along lines of constant entropy per baryon. 

We remark that most of the work done in this area has been performed with non strange 
quarks only and, when strange quarks are considered, mixing between the flavors u, d, and 
s has not been taken into account Our SU(3) version of the NJL model includes a term 
that incorporates the axial anomaly of QCD, and is responsible for the mechanism of flavor 
mixing. We relate the discontinuity of the order parameter to other discontinuities of physical 
quantities such as, for instance, the entropy. We are particularly interested in confronting 
our calculation, in what concerns to the notion of a second order phase transition due to 
nonvanishing current quark masses, with those of any classical mean field theory. From 
lattice calculations it is well known that the strange quark mass plays a decisive role in the 
location of the CEP. 

On the other hand, information on the nature of excitations and the strength of their 
interaction in the QGP would be crucial in the experimental search. Also in this context it 
is relevant to confront first-principle based approaches with the results of phenomenological 
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models like the NJL model. 

We organize the work in four main steps. First, after the presentation of the model 
formalism (Sec. II), we discuss the behavior of the equations of state and analyze the chiral 
phase transition (Sec. III). The well known universality hypothesis of phase transitions will 
be considered. Second, we study the behavior of relevant physical quantities in the T — /i^ 
plane (Sec. IV). Third, we analyze the phase diagrams in the T — fiB plane looking for the 
location of the critical end point and the behavior of susceptibilities (Sec. V). Finally, we 
discuss signs of partial and effective restoration of chiral symmetry (Sec. VI), looking for 
the convergence of chiral partners. We conclude in Sec. VII with a brief summary of our 
results. 



II. FORMULATION OF THE MODEL 



The Lagrangian of the SU(3) NJL model 



21 



22l | is given by 



a=0 



+ 9D 



det [q{l + 75)g] + det [q{l - 75)5] 



(1) 



The column vector q = {u, d, s) represents the quark field with three flavors, Nf = 3, and 



three colors, Nc = 3. A" are the Gell-Mann matrices, a = 0, 1, . . . , 8, A 
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The Lagrangian ([T]) is invariant under chiral SUL(3)®SUi?(3) transformations if we put 
rrii = 0, where rrii are the current quark masses (m = diag(m„, m^, m^)). The last term in 
([T]) breaks the Ua(1) symmetry. This term is a reflection of the axial anomaly in QCD. 

The model Lagrangian ([T]) can be put in a form suitable for the bosonization procedure 
after an adequate treatment of the last term, allowing to obtain a four quark interaction 
from the six quark interaction. Then the following effective quark Lagrangian is obtained: 

C,jf = q{ii^d, - m)q + ( g A" g )(g )] + P,5[( g ^ 75 A" g ) ( gz 75 A^ ) ], (2) 

where the projectors Sab , Pab are given by: 



Sab = gs^ab + gDP>abc {q^q) , 

Pab = 9sSab - QoDabc {q>^''q) ■ 



(3) 
(4) 
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The constants Dabc coincide with the SU(3) structure constants dabc for a, 6, c = (1,2,..., 8) 
and Doab = ~^^ab, -Dooo = -y/f- The hadronization procedure can be done by the integra- 
tion over the quark fields in the functional integral with ([2]). The natural degrees of freedom 
of low-energy QCD in the mesonic sector are achieved which gives the following effective 
action: 



iTi In 



trd^-m + aa\'' + it^5){^aXn . (5) 



The notation Tr stands for the trace operation over discrete indices {Nf and Nc) and inte- 
gration over momentum. The fields a"" and ip°' are scalar and pseudoscalar meson nonets, 
respectively. 

The first variation of the action (JSj) leads to the gap equations. 

Mi = nii - 2gg{qiqi) - 2g^{qjqj){qkqk) , (6) 

with k = u,d,s cyclic. Mj are the constituent quark masses and the quark condensates 
are given by: {(jiqi) = — iTr[5'i(p)], Si{p) being the quark Green function. 

The baryonic thermodynamic potential of the grand canonical ensemble, Q{T,V, fii), is 
also obtained directly from the effective action (I5l). So we take the temperature T, the volume 
V and the chemical potential of the z-quark (yUj) as the full independent state variables. 

The relevant equations of state for the entropy S, the pressure p, and the particle number 
Ni, as well as the internal energy E, follow from well known expressions like the Gibbs- 
Duhem relation 

Q{T, V,n,) = E-TS- J2 f^^^ ■ (7) 

i=u,d,s 

The following expressions are obtained: 

E = - — V 2^ < P dp — (1 - Tii - Hi) I 

i=u,d,s ^ * ^ 

~9s ^ {{qiqi})'^ -2gD{uu){dd){ss) , (8) 



5" = --^V J P^dp I Klnnj + (1 - rii) ln(l - rii)] + [rii ^ 1 - ni\ \ , 



(9) 
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Ni = :^v I p^dp im - n,) . (10) 

V is the volume of the system and the quark density is determined by the relation pi = Ni/V. 
In the previous equations n.j and Ui are the quark and antiquark occupation numbers 

1 1 



Hi 



(11) 



We define Pb = |(At« + /Wd + /Ws) and the baryonic matter density as ps = |(Pm + Pd + Ps)- 
As usual, the pressure and the energy density are defined such that their values are zero in 
the vacuum state [23|: 

1 ^^2^ 



V 



[Q{p„T)-Q{0,0)] 



e{p^,T) = -[E{p,,T)-E{0,0)]. 



(13) 



The baryon number susceptibility is the response of the baryon number density Pb{T, pi) 
to an infinitesimal variation of the quark chemical potential pi 2J] : 



Xb 



dpi 

dpiy J, 



(14) 



Another relevant observable, in the context of possible signatures for chiral symmetry 
restoration in the hadron-quark transition and in transition from hadronic matter to the 



QGP 
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261], is the specific heat which is defined by [jj 



C 



V [df 



T 
V 



dS_ 
df 



{dNi/dp. 



i T 



(15) 



where we have transformed the derivative {dS/dT)^^ using the formula of the Jacobian. 
In fact, we work in the grand canonical ensemble where {T,V,pi) are the set of natural 
independent variables (still holding Ni and V fixed). 

By expanding the effective action (l5l) over meson fields, we get an effective meson action 
from which we can obtain the meson propagators. In the present work we are only concerned 
with 7r° and a mesons. Starting with the pseudoscalar mesons we have the effective meson 
action: 



where n^j(P) is the polarization operator. 



nr,(P) = iN^j ^trz, [5,(p)(A'^),,(z75)5,(p + P)(A%(z75)] 



(16) 



(17) 



with tTD is the trace over Dirac matrices. The expression in square brackets in (fT6!) is the 
inverse non-normahzed meson propagator (D^(P))~^. 
The inverse meson propagator for is given by 

D;,}iP) = i-P^oj^,iP), (18) 

with 

Ptto = gs + 9D {qsQs) (19) 

and where the polarization operator of the 7r° meson takes the form 

JuuiPo) = 4 [2/r - Jr(Po)] . (20) 

The integrals !{ and l2{Po) are given in Appendix A. 

The mass of the 7r° meson can be determined by the condition D^o{M^o, 0) = and the 
quark-meson coupling constant is evaluated as 

The procedure to describe scalar mesons is analogous. We present below the most relevant 
steps. Keeping now the scalar mesons only in ([5]), we have the effective meson action 

<}M = -l^" K - nf.(P)] = -^a«(I^f,(P))^V^ (22) 

with nf^(P) being the polarization operator, which in the momentum space has the form of 
( TT7|) with (275) substituted by the identity matrix. 

To consider the a meson we take into account the matrix structure of the propagator in 
f l22|) . For the isospin symmetry considered in the present work {quQu) = {QdQd), and the 
matrices Sab and Ilf^ are reduced to 

S.-f^" ° ] and ni^f"- ° V (23) 

where the matrix elements are given in Appendix A. 

The mass of the a meson can be determined by the condition D~^{M„, 0) = 0, where 

D;^ = {A + C)- ^{C - Ay + (24) 

with A = Ss8- Anfo(P), C = Soo- Anf8(P), B = -{S08 + AUl^iP)) and A = SooSss - 



Finally, the model is fixed by the coupling constants gs and gr,^ the cutoff in three- 
momentum space A, which is used to regularize the momentum space integrals and the 
current quark masses rrtj. For numerical calculations in physical conditions we use the 



parameter set 
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27 



28 



29j: 



TUd = 5.5 MeV, m, = 140.7 MeV, gsA'^ = 3.67, 



— 12.36 and A = 602.3 MeV, that has been determined by fixing the values M^r = 135.0 
MeV, Mk = 497.7 MeV, = 92.4 MeV, and M^/ = 960.8 MeV. For the quark condensates 
we obtain: (g„ g„) = {qaQd) = -(241.9 MeV)^ and = -(257.7 MeV)^ and for the 

constituent quark masses M„ = = 367.7 MeV and Mg = 549.5 MeV. 



III. EQUATIONS OF STATE AND PHASE TRANSITION 

We will start the discussion of the phase diagram of the NJL model ([1]) by analyzing the 
behavior of the pressure/energy per particle as a function of the baryonic density, paying 
special attention to the Gibbs criteria. Our model of strong interacting matter can simulate 
either a region in the interior of a neutron star or a hot and dense fireball created in a 
heavy-ion collision. In the present work we focus our attention in the last type of systems, 
so we impose the condition /ig = 0, since electrons and positrons are not involved in the 
strong interaction. So, we naturally get the chemical equilibrium condition fiu = f^d = f^s = 
fiB that will be used. This choice allows for equal constituent quark masses = 
and approximates the physical conditions at RHIC. In this respect, we remind that in a 
relativistic heavy-ion collision of duration of ~ 10~^^ s, thermal equilibration is possible only 
for processes mediated by the strong interaction rather than the full electroweak equilibrium. 

Let us discuss our results for the pressure/energy per baryon at zero temperature that 
are plotted in Fig. [T]as a function of Pb/ Po (solid lines), where po = 0-17 fm~'^ is the normal 
nuclear matter density. The pressure has three zeros, respectively, at ps = 0, 0.43po) 2.36po, 
that correspond to the extreme of the energy per particle. For Pb < 0.2po the pressure 
and compressibility are positive, so the system can exist in a uniform gas phase, but it will 
not survive indefinitely, since the zero density state is energetically favored; for 0.2po < 
Pb < 0.43po the system is unstable since the compressibility is negative, in fact Pb = 0.43po 
corresponds to a maximum of the energy per particle; for 0.43po < Pb < 2.36po; the pressure 
is negative, and the third zero of the pressure, pb = 2.36po, corresponds to an absolute 
minimum of the energy (see Fig. [T] (right panel)). The appearance of an absolute minimum 
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FIG. 1: Pressure (left) and energy per particle (right) as a function of the density at different 
temperatures. The points A and B (left panel) illustrate the Gibbs criteria. Only in the T = 
line the zero-pressure point is located at the minimum of the energy per particle. 



of the energy means the possibihty for finite droplets to be in mechanical equilibrium with 
the vacuum at zero -pressure [P = 0). Above pb = 2.36po; which we define as p^, we have 
again a uniform gas phase. So, for densities < ps < Pb the equilibrium configuration 
is a mixed phase. This is because the Gibbs criterion of equal P and /i_B is satisfied and, 
therefore, the phase transition is a first order one: the thermodynamic potential has two 
degenerate minima at which two phases have equal pressure and chemical potential and 
can coexist. Such a situation is possible in regions where the gap equations have several 
solutions for the quark masses. 

Summarizing the results at T = 0, the behavior described allows the following interpreta- 
tion: the uniform nonzero density phase will break up into stable droplets with zero pressure 
and density = 2.36po in which chiral symmetry is partially restored, surrounded by a 
nontrivial vacuum with ps = P = (see also 10|, l23|, |27|, |30|, l3l(|)- In fact, for our choice of 



23 



where 



the parameters the critical point at T = satisfies to the condition p^ < 
M^"''^ is the mass of the i-quark in the vacuum. This can be seen by comparing p^' = 361 
MeV (see the T-axis of Fig. 2, left panel) with the quark masses M^"-'^ = M^"-'^ = 367.7 
MeV and MJ'^'= = 549.5 MeV. 

As can be seen from Fig. [H as the temperature increases, the first order transition persists 
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FIG. 2: Phase diagram in the SU(3) NJL model. The left (right) part corresponds to the T — fis 
{T — pb) plane. Solid (dashed) line shows the location of the first order (crossover) transition. The 
dotted lines shows the location of the spinodal boundaries of the two phase transitions (shown by 
shading in the right plot). 

up to the CEP. At the CEP the chiral transition becomes of second order. Along the hne of a 
first order phase transition the thermodynamic potential has two degenerate minima. These 
minima are separated by a finite potential barrier making the potential nonconvex. The 
height of the barrier is largest at zero temperature and finite quark chemical potential and 
decreases towards higher temperature. At the CEP the barrier disappears and the potential 
flattens. This pattern is characteristic of a first order phase transition: the two minima 
correspond, respectively, to the phases of broken and restored symmetry. The borders of 
the coexistence marked by the dotted lines in Fig. [2l The domain between the 

two dotted lines has metastable states which are characterized by large fluctuations. They 
are also solutions of the gap equations but their thermodynamic potential is higher than 
for the stable solutions. The left dotted curves represent the beginning of the metastable 
solutions of restored symmetry in the phase of broken symmetry, while the right dotted 
curves represent the end of the metastable solutions of broken symmetry in the restored 
symmetric phase. We also represent in Fig. [2] (right panel) the region where the solutions 
of the gap equations are unstable. 

The location of the CEP is found to be at T^-^^ = 67.7 MeV and pg^^ = 1.68po 
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(fJ'B^^ = 318.4 MeV). For temperatures above the CEP the thermodynamic potential has 
only one minimum and the transition is washed out: a smooth crossover takes place. 

Finally, we will focus again on the energy per baryon. In Fig. [T] (right panel), we plot 
the density dependence of the energy per baryon at different temperatures. We observe that 
the two points, zero of the pressure and minimum of the energy density, are not the same 
at finite temperature. In fact, as can be seen from Fig. [1] (left panel), states with zero 
pressure are only possible up to the maximal temperature ~ 38 MeV. We notice that 
the zero-pressure states persist up to temperatures of 70 MeV in a two-flavor NJL model 
where equal chemical potentials of quarks and antiquarks is assumed 30|]. For T < Tm the 
zero-pressure states are in the metastable density region and, as soon as T 7^ 0, they do not 
coincide with the minimum of the energy per particle. 

The arguments just presented allow to clarify the difference between confined quark 
matter (in hadrons) and bounded quark matter (droplets of quarks). As would be expected, 
the binding mechanism is weaker than the confining one (nonexistent in the NJL model). 
As a matter of fact, in spite of the existence of a binding energy for the droplets of quarks 
at T = 0, we verify that it is not possible to avoid the evaporation of the bounded quarks 
for arbitrarily small temperatures. 

More detailed information concerning the structure of the phase diagram will be given in 
Sec. V. 



IV. THERMODYNAMIC QUANTITIES IN THE T - hb PLANE 

For a better understanding of the thermodynamics of the phase transitions, we analyze 
in this section the behavior of the thermodynamical quantities that are the most relevant 
to discuss the physics across the first order phase transition. With these quantities, we can 
also discuss the latent heat which is inherent to this phase transition. 

The pressure is plotted in the left- hand side of Fig. [3] (upper part), which shows a 
continuous behavior for all points of the phase diagram. In a first order phase transition a 
discontinuity occurs in the first derivatives of the pressure (or the thermodynamic potential) 
with respect to and T, i.e., the baryon number density and the entropy density, respec- 
tively. In fact, as can be seen in the right side of Fig. [31 the entropy density is discontinuous 
in the first order phase transition region (T < T*^^^, fis > fJ'B^^)- A similar behavior is 
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FIG. 3: Pressure (left side of upper part), entropy density (right side of upper part) and energy 
density (down part) as functions of the temperature and the baryonic chemical potential. 

found for the energy density, whose curves show that the first order phase transition, strong 
at T = 0, turns into a less abrupt one as the temperature increases (see Fig. [3|, lower part). 
In the crossover transition (T > T'^^^, fiB < f^%^^) the thermodynamic quantities change 
rapidly within a narrow range of values of T and /i^, but the pressure and all its derivatives 
remain continuous, as shown in Fig. [3l 

The discontinuities of the entropy and energy densities disappear at the CEP, which 
location can not be determined by universality arguments. The same is not true concerning 
local singular behavior of thermodynamic quantities around the CEP that will be discussed 
in the next section through the critical exponents. 

Let us now analyze what information concerning the latent heat we can get from our 
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results. As already referred in Sec. Ill, along the line of first order phase transition the 
thermodynamic potential has two degenerate minima that are separated by a finite barrier. 
This barrier is largest at zero temperature and finite chemical potential and decreases to- 
wards higher temperature. At the CEP the barrier disappears, which means that there is 
no latent heat at this point. 

As a grand canonical approach is applied to our model of strong interacting matter, 
the independent quantities T and fiB represent the state variables which can be externally 
controlled. So, the conjugate of the intensive variables T and /is in the Legendre trans- 
formation — the entropy density s and the baryonic density ps — provide a more natural 
description. By analyzing first the gap in the curves of the entropy (Fig. [31 right side of 
upper part), we see that the latent heat decreases for small temperatures, which is not the 
expected behavior. This analysis is, however, not sufficient; both the baryonic density and 
the entropy density contributions should be examined for more reliable information about 
the latent heat. We remember that the gap of the baryonic density across the first order 
phase transition is largest at zero temperature and finite chemical potential and vanishes at 
the CEP (see Fig. El right panel). The discontinuities in the energy density include both the 
entropy and the baryonic density contributions and, as can be seen in Fig. |3l they display 
the expected behavior: the latent heat increases for decreasing temperatures. 

Finally, to understand the thermodynamics of matter created in relativistic heavy-ion 
collisions, it is convenient to calculate thermodynamic quantities along lines of constant 
entropy per baryon number, the so-called isentropic lines. Most of these studies have been 



done on lattice calculations for two-flavor QCD at finite fiB 33|], where nonphysical mass 



spectrum that corresponds to a too large of pion mass ~ 770 MeV, has been used. Such 
studies predict that the effects of the CEP change only slowly as the collision energy is 
changed as a consequence of the attractor character of the CEP (l2 |. 

Our model calculations for the isentropic lines in the T — fis plane are shown in Fi g. [4l 



The behavior we find is somewhat different from those claimed by other authors 33|, l34, l35| , 
where a phenomena of focusing of trajectories towards the CEP is observed. 

The isentropic trajectories in the phase diagram (Fig. |4]) indicate that the slope of the 
trajectories goes to large values for large T. This behavior is related to the effects of the 
strange quark mass in our model. In fact, at high temperatures the relation ps > Mg is 
verified, allowing for a more pronounced decrease of Ms 271]. Although the entropy and the 
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FIG. 4: Two perspectives of the entropy per baryon number in the T — fiB plane. The vicinity of 
the CEP is enlarged in the right panel. 

baryon number density, at high temperatures, are sensitive to the regularization procedure 
used jsol . Is?! , this effect is not relevant for the present situation. The same is not true with 
respect to the effects of the value of the cutoff itself in the regime of low temperatures as 
will be shown below. 

In a small range of s/pb around 0.7 (see Fig. |U right panel), we observe a tendency of 
convergence of these isentropic lines towards the CEP. These lines come from the region of 
symmetry partially restored in the direction of the crossover line. For smaller values of s/ps, 
the isentropic lines turn about the CEP and then attain the first order transition line. For 
larger values of s/ps the isentropic trajectories approach the CEP by the region where the 
chiral symmetry is still broken, and also attain the first order transition line after bending 



toward the critical point. As already pointed out in [32], this is a natural result in these 
type of quark models with no change in the number of degrees of freedom of the system 
in the two phases. As the temperature decreases a first order phase transition occurs, the 
latent heat increases and the formation of the mixed phase is thermodynamically favored. 

Finally, we remark that all isentropic trajectories directly terminate at T = at the first 
order transition line, without reheating in the mixed phase as verified in the "zigzag" shape 



of 



33 



M 



35 



It is also interesting to point out that, in the limit T — > 0, it is verified that 



s — ^- and ^ 0, as it should be. This behavior is in contrast to 321] (right panel Fig. 9) 
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using the NJL model in the SU(2) sector and is related to our more convenient choice of the 
model parameters, mainly a lower value of the cutoff. This can be explained by the presence 
of droplets at T = whose stability is quite sensitive to the choice of the model parameters. 
In fact, as referred in Sec. Ill, our choice of the parameters has important effects: we verify 
that, at T = 0, the phase transition connects the vacuum state (P = 0,pb = 0) directly 
with the phase of chiral symmetry partially restored (P = 0, p = p^) and the critical point 
of the phase transition in these conditions satisfies to fii < M^""^, where M^"-'^ is the mass 



of the i-quark {i = u, d, s 



of the quark droplets 
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in the vacuum. This condition fulfills the criterium of stability 



271]. In addition, it is also crucial to the satisfaction of the third 



law of thermodynamics in the limit T ^ 0. This cutoff effect has an identical role in the 
formation of stable droplets on both SU(2) and SU(3) NJL models. 



V. PHASE DIAGRAMS AND SUSCEPTIBILITIES IN THE VICINITY OF THE 
CRITICAL END POINT 

In this section we analyze with more detail the phase diagrams in different conditions in 
the T — fiB plane. Lattice QCD calculations have established the transition to a phase where 
quarks and gluons are deconfined at temperatures larger than ~ 150 MeV and zero baryon 
density. Depending on the number of quark flavors Nf = 2 or Nf = 3, and on the masses of 
the quarks, different situations can occur and the transition from hadronic matter to QGP 
may be of first order, second order, or a crossover transition. To confront the model results 
with the universality arguments, we will discuss the class of the critical points by changing 
the current quark masses in SU(2) and SU(3) versions of the NJL model. 



A. Characteristics of the T — phase diagram 

We start by analyzing the differences between the three-flavor NJL model and its simpler 
version in the SU(2) sector. The phase diagrams for both models are presented in Fig. Oas 
a function of fis and T. 

Concerning the SU(2) model, and using physical values of the quark masses: m„ = = 
5.5 MeV, we find that the CEP is localized at T^'^^ = 79.9 MeV and /ig^^ = 331.72 MeV 
(Pb^^ = 2.26po)- We also verified that, in the chiral limit, the transition is of second order 
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FIG. 5: Phase diagram in the SU(2) (left) and SU(3) (right) NJL models. The solid line represents 
the first order phase transition, the dashed line the second order phase transition, and the dotted 
line the crossover transition. The size of the critical region is also plotted for several values of 
Xb/Xb^'^ ■ '^^^ TCP in the right panel is found for ?n„ = = MeV and vtls = 140.7 MeV. 

at fiB = and, as fiB increases, the line of second order phase transition will end in a first 
order line at the TCP. The TCP is located at //^^^ = 286.1 MeV and T^^^ = 112.1 MeV. 

For the SU(3) NJL model, also in the chiral limit (m„ = rrid = nia = 0), we verify that 
the phase diagram does not exhibit a TCP: chiral symmetry is restored via a first order 
transition for all baryonic chemical potentials and temperatures (see right panel of Fig. [5]). 
According to lattice analysis, this pattern of chiral symmetry restoration should remain even 
when the strange quark acquires a nonzero current mass, provided it is lower than a critical 
value {irts < m^"*), and niu = rrid = is still kept. The value for m^"* is not settled yet. 



those found in lattice 39|] or in model calculations 40|, |41| being lower than the physical 



strange current quark mass (m^ ^ 150 MeV). We found m^"* = 18.3 MeV in our model 15| . 



ower than lattice values 



4l|. 



391] but consistent with what it is expected in these type of models 



When nis > rrig"^, at /i^ = 0, the transition is of second order and, as /i^ increases, the 
line of the second order phase transition will end in a first order line at the TCP. The TCP 



for nis = 140.7 MeV is the closest to the CEP [15j and is located at = 265.9 MeV 



16 



and T^'^^ = 100.5 MeV. If we choose = 7^ 0, instead of second order transition we 



have a smooth crossover whose critical 
physical values for the quark masses 



line will end in the first order line at the CEP. Using 
28j: mu = md = 5.5 MeV, = 140.7 MeV, this 
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point is localized at T^^^ = 67.7 MeV and /ig^^ = 318.4 MeV (pg^^ = 1.68po)- 

We point out that both situations are in agreement with what is expected at /i^ = 0: the 
chiral phase transition at the chiral limit is of second order for Nf = 2 and first order for 



Nf>3 



We also observe that the critical region is heavily stretched in the direction of the crossover 
transition line, in both Nf = 2 and Nf = 3 shown in Fig. [5l To estimate the 

critical region around the CEP we calculate the dimensionless ratio Xb/x^b'^^ ^ where x^b'^'^ is 
the chiral susceptibility of a free massless quark gas. The left (right) panel of Fig. O shows a 
contour plot for two fixed ratios Xb/x'b'^'^ ~ 2-0; 3.0 in the phase diagram around the CEP. 



B. Behavior of xb and C in the vicinity of the critical end point and their critical 
exponents 

The phenomenological relevance of fluctuations in the finite temperature and chemical 
potential around the CEP/TCP of QCD has been recognized by several authors. If the 
critical region of the CEP is small, it is expected that most of the fluctuations associated 
with the CEP will come from the mean field region around the CEP [l^. The size of the 
critical region around the CEP can be found by calculating the baryon number susceptibility, 
the specific heat and their critical behaviors. 

To a better understanding of the critical behavior of the system, we also analyze in some 
detail what happens in the SU(2) case, sector to which there is more information in the 
literature 4^. 



As is well known, the 
at T = T^^^ 
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Daryon number susceptibility, xb, and the specific heat, C, diverge 
44| . In order to make this statement more precise, we will focus on 
the values of a set of indices, the so-called critical exponents, which describe the behavior 
near the critical point of various quantities of interest (in our case e and a are the critical 
exponents of Xb and C, respectively). The motivation for this study arises from fundamental 
phase transition considerations, and thus transcends any particular system. These critical 
exponents will be determined by finding two directions, temperature and magnetic-field- 



17 



Quantity critical exponents/path SU(2) NJL 



SU(3) NJL Universality 



e/ ^« 0.66 ±0.01 0.67 ±0.01 2/3 

XB e' / 0.66 ±0.01 0.68 ±0.01 2/3 

7ij/^. 0.51 ±0.01 0.49 ±0.02 1/2 

• a = 0.59 ±0.01 0.61 ±0.01 2/3 

a I 

t ai= 0.45 ±0.01 — — 

i 

C oil 0.69 ±0.01 0.67 ±0.01 2/3 



a I 0.40 ±0.01 0.45 ±0.02 1/2 

T 



TABLE I: The arrow ^ • ( | indicates the path in the [Ib {T)— direction to the CEP (TCP) 

T 



for fiB < /ig^^ {T < T^^^). 



like, in the T — hb plane near the CEP, because, as pointed out in |45|, the strength 
of the divergence is governed by the critical exponents whose values depend on the path 
approaching the CEP. 

Considering the baryon number susceptibility, if the path chosen is asymptotically parallel 
to the first order transition line at the CEP, the divergence of xb scales with an exponent 
7^. In the mean field approximation it is expected that 75 = 1 for this path. For directions 
not parallel to the tangent line the divergence scales as e = 2/3. These values are responsible 
for the elongation of the critical region, xbi being enhanced in the direction parallel to the 
first order transition line (see Fig. E]). 

To study the critical exponents for the baryon number susceptibility (Eq. \\M we will 
start with a path parallel to the /xs-axis in the T — fiB plane, from lower fiB towards the 
critical at fixed temperature T = T^^p _ Using a linear logarithmic fit 

InxB = -eln l/iB - yug-^^l + ci, (25) 

where the term ci is independent of /i^, we obtain e = 0.67 ± 0.01, which is consistent with 
the mean field theory prediction e = 2/3. 

We also study the baryon number susceptibility from higher [ib towards the critical 
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h'b^^ . The logarithmic fit used now is Inx^ = — e'ln|/iB — + c'^. Our result shows 

that e' = 0.68 ±0.01 which is very near the value of e. This means that the size of the region 
we observe is approximately the same independently of the direction we choose for the path 
parallel to the /i^-axis. These critical exponents, calculated in both SU(2) and SU(3) NJL 
models, are presented in Table 1. 

For comparison purposes with the universality /mean field predictions, the calculated 
critical exponents at the TCP are also presented in Table I. It is found that the critical 
exponent for xb, 1b once we are in the TCP, has the value 7^ = 0.49 ± 0.02, for the SU(3) 
NJL model and 7^ = 0.51 ± 0.01, for the SU(2) NJL model. These results are in agreement 
with the mean field value (7^ = 1/2), and show that the behavior of the baryon number 
susceptibility is similar in both SU(2) and SU(3) versions of the model. 

Paying now attention to the specific heat (Eq. [T5|) around the CEP, we have used a path 
parallel to the T-axis in the T — fiB plane from lower/higher T towards the critical T^^^ 
at fixed = A^b^^- In Fig- E] (upper part) we plot C as a function of T close to the 
CEP in a logarithmic scale for both SU(2) and SU(3) calculations. In this case we use the 
linear logarithmic fit InC = —a In \T — T^^p\ -|- where the term C2 is independent of T. 
Starting with the SU(2) case, we observe in the left panel that, for the region T < T'"^^, 
we have a slope of data points that changes for values of |T — T'^'^-^l around 0.3 MeV. We 
have fitted the data for |T - T^^^| < 0.3 MeV and |T - T^^^| > 0.3 MeV separately and 
obtained, respectively, the critical exponent a = 0.59 ± 0.01 and ai = 0.45 ± 0.01, which 
lave a linear behavior for several orders of magnitude (see also Table I). As pointed out in 
isl ]. this change of the exponent can be interpreted as a crossover of different universality 
classes, with the CEP being affected by the TCP. It seems that in our model the effect of 
the hidden TCP on the CEP is relevant for the specific heat contrarily to what happens to 
Xb- 

We also observe that there is no clear evidence of the change of the slope of the fitting of 
data points in the three- flavor NJL model (see Fig. [6], right panel of upper part, and Table 
I). In fact, now we only obtain a critical exponent a = 0.61 ± 0.01 when the critical point 
is approached from below. When the critical point is approached from above the trivial 
exponent a' = 0.67 ± 0.01 is obtained. 



To explore the possible effect of the hidden TCP on the CEP, as suggested in Refs. 



we analyze the behavior of the specific heat around the TCP. As shown in Fig. O (lower part) 
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FIG. 6: Upper part: Specific heat as a function of \T — T^^P\ at the fixed chemical potential 
^CBP £qj, gu(^2) (left) and SU(3) (right) NJL models. Lower part: Specific heat as a function of 
\T - T^CP^ ^j^g g^g^i chemical potential fi^^^ for SU(2) (left) and SU(3) (right) NJL models. 

and Table I, we find nontrivial critical exponents a = 0.40 ± 0.01 and a = 0.45 ± 0.01, for 
SU(2) and SU(3) cases, respectively. This result, in spite of being close, is not in agreement 
with the respective mean field value (a = 1/2). However, they can justify the crossing effect 
observed. We notice that the closest distance between the TCP and the CEP in the phase 
diagram occurs in the T-direction ((T"^*"^ — T*"^^) < — jj^^)). 
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The inconsistency with the mean field values only occurs for the exponent a as can be 
seen from Table I. According to what was suggested by universality arguments in [l^ , it was 
expected that xb and C should be essentially the same near the TCP and the CEP, which 
would imply a = e = 2/3at the CEP. Nevertheless we observe that the nontrivial values 
of a in the TCP and in the CEP are consistent within the NJL model for both, SU(2) and 
SU(3) versions of the model, and they refiect the effect of the TCP on the CEP. We also 
stress that the universality arguments are so general that they give no quantitative results 
and, due to the lack of information from the lattice simulations, they should be confronted 
with model calculations. Our results seem particularly interesting because the NJL model 
shares with QCD some features, such as the dynamics of chiral symmetry. In particular, 
the physics underlying the critical singularities in the QCD diagram is associated with this 
fundamental property of strong interaction. So, the NJL model is an useful framework 
allowing for insights to the difficult task of the analysis of the QCD phase diagram at finite 
temperature and chemical potential. 

The eventual difference between the values of the C and xb critical exponents can be 
interesting in heavy-ion collisions experiments. 



VI. PARTIAL AND EFFECTIVE RESTORATION OF CHIRAL SYMMETRY 



As we have shown in previous sections, thermodynamics provides a well- established 
procedure, as for instance the Gibbs criterion, to determine the critical points for the phase 
transition in the first order region. It follows that these critical points are signalized by the 
discontinuity of several relevant observables (masses, quark condensates) at some critical 
chemical potential, a situation that does not happen in the crossover region, where these 
observables are continuous. At present, the criterion most commonly accepted, and that 
will be used here, to define the critical point in the crossover region, is to identify this point 
as the infiection point of the quark masses d'^M/dT'^ = j23| or, equivalently, of the quark 



condensates 
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47|. This criterion is numerically equivalent to the one first proposed by M. 



Asakawa and K. Yazaki that defines the point where the constituent quark masses decrease 
to half of their values in the vacuum (M^ = Mu(0)/2) 48|], as the critical point. From this 
point on the quark masses decrease quickly. 

Both in the first order and in the crossover regions it is verified that the quark masses. 
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especially for the non strange quarks, decrease strongly at the critical point. However, at 
this point different observables violating chiral symmetry are still far from zero, like the 
quark condensates, the pion decay constant, and the difference between the masses of the 
chiral partners. One can say, therefore, that at the critical point there occurs only a partial 
restoration of chiral symmetry. 

In view of what was said above we use the following criteria: we define the point in the 
T — fiB plane for the phase transition associated with partial restoration of chiral symmetry 
as the inflexion (discontinuity) point for the quark masses, and define the point for effective 
restoration of chiral symmetry as that one where the masses of chiral partners become 
degenerate. This is also signaled by the merging of the 7r° and a spectral functions 49 1. 

As we include the strange sector in this study, the consequences of the nonvanishing 
anomaly term (mixing effects) on the strangeness content of mesons and mixing angles must 
be analyzed. In fact, as the temperature (density) increases, the mixing ang les get close to 



their ideal values and the strangeness content of the mesons change [28|, l29|] : the masses of 
the mesons that become almost non-strange, a and rj, converge, respectively, with those of 
the non strange mesons 7r° and Oq, while that of the rj', that becomes essentially strange, 
does not get close to /o (see 29|] (Fig. 2, Case I)); the convergence of the chiral partner 
{k,K), that has a us structure, occurs at higher temperatures and is probably slowed by 
the small decrease of the constituent strange quark mass, Mg. For the purpose of discussing 
the effective restoration of chiral symmetry, we restrict our analysis to the chiral partners 
(vr^, 0") that behave in a qualitatively similar manner as the pair (ao,?7). 

The behavior of the masses of the chiral partners (vr^, a) in the limiting cases (T 7^ 0, 
Pb = 0) and (T = 0, ps ^ 0) are qualitatively similar and well known from the literature: 
they both converge at a certain value of the temperature (density). The main difference 
between the finite temperature and the finite density case is that, in the first one, the 
degeneracy of the chiral partners occurs in a range of temperatures where the mesons are 
no longer bound states: the 7r° dissociates in qq pair at the Mott temperature Tm-kO = 212 



MeV 
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271], and the a at the Mott temperature T/v/o- = 160 MeV; for the finite density 



case, the mesons are always bound states. 

Interesting information can be obtained by calculating the masses of the 7r° and a mesons 
as a function of T and Pb{pb) which allows us to obtain a curve in the T — pb{pb) plane. 
This curve defines the line where the mesons became degenerate (Fig. [7]). In Fig. [7] we 
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FIG. 7: The effective restoration of chiral symmetry, the phase transition and the Mott hnes for 
■K^ and a mesons in the T — p_b(/U_b) plane. 

also represent the "Mott lines" for the ti^ and the a, as well as the critical line. As can 
be seen, the phase transition associated to partial restoration of chiral symmetry occurs 
above the Mott line for the pion and below the Mott line for the sigma, in most of the first 
order phase transition region; the opposite happens in the crossover region. Concerning 
the effective restoration of chiral symmetry, one can see, from the line of convergence of 
the chiral partners, that it happens after the partial restoration of chiral symmetry and the 
dissociation of the two mesons. 

As we already saw, there are dramatic changes in the behavior of some thermodynamic 
functions such as the specific heat and the quark number susceptibilities around the CEP. 
So, due to their role as signals for restoration of chiral symmetry it is demanding to discuss 
the behavior of the chiral partners (tt^, a). 

First we notice, in Fig. [TJ that the two Mott lines cross in the first order region at a point 
just bellow the CEP. This is probably a remnant of the situation in the chiral limit where 
the transition is second order and the pion and sigma dissociate at the same point. 

In Fig. [8] we plot the pion and sigma masses as functions of the baryonic chemical 
potential for three different temperatures: T = 40 MeV < T^^^^^ j^cep _ Q^j rj j^^y 
T = 100 MeV > T^^P. For T = 40 MeV and jje ~ 350 MeV, a discontinuity is visible 
in the evolution of the masses, signaling a first order phase transition. However, according 
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FIG. 8: Degeneracy of the chiral partners (vr", a) for different temperatures around the CEP. 

to our criterion, the effective restoration of the chiral symmetry only happens at /i^ ~ 380 
MeV. At the CEP (T = 67.7 MeV; hb = 318.4 MeV), the sharp decrease (increase) of the 
sigma (pion) meson masses reflect the nature of the second order phase transition. Once 
again the effective restoration of chiral symmetry only happens at ^ 370 MeV. When 
T = 100 MeV > T*^^^, we have a crossover and the meson masses have a smooth behavior. 
In this case, the effective restoration of the chiral symmetry happens at /x^ ~ 355 MeV. 



VII. CONCLUSIONS 



The properties of the QCD transition at vanishing chemical potential depend on the 
number of quark flavors and on their masses. The critical temperatures of Tc ~ 155 MeV or 
as high as Tc ~ 260 MeV have been reported in the literature. Presently, considering also 
nonvanishing chemical potential, some lattice calculations locate the CEP a.t T ^ 160 MeV 
and /i_B ~ 360 MeV 6|. However, the existence and location of the CEP are not conclusive 



even for lattice calculations 



50|. 



We proved that our model calculation has been able to reproduce the qualitative phase 
structure features, and we also obtain the location of the CEP. We have obtained, at zero 
baryon chemical potential in the SU(3) NJL model, values for the critical temperature 
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around 120 — 200 MeV. The transition is first order in the chiral hmit (m„ = rrid = = 0). 
Furthermore, when m„ = = and > mf'** {mf^^ = 18.3 MeV) the transition is of 
second order ending in a first order hne at the TCP. Finally, when also = rrid ^ 0, there is 
a crossover for all values of and the location of the CEP depends strongly on the strange 
quark mass. Contrarily to what happens in the three-flavor NJL model, we find a TCP in 
the two-flavor NJL model in the chiral limit. This agrees with what is expected at fiB = 0: 
for rrii = the chiral restoration happens via a second order phase transition for Nf = 2, 
and via a first order for A'^^ > 3. For realistic values of the current quark masses the CEP 
is located at T^^^ = 79.9 MeV and /ig^^ = 331.7 MeV for Nj = 2, and at T^^^ = 67.7 
MeV and /ig^^ = 318.4 MeV for Nf = 3. 

The pattern characteristic of a first order phase transition has also been analyzed through 
several equations of state and the latent heat. For example, we verified that states (droplets) 
in mechanical equilibrium with the vacuum state at P = are found at zero temperature. 
This leads to nontrivial consequences for the behavior of the isentropic trajectories which 
terminate at T = at the first order transition line. Our convenient choice of the model 
parameters, which allows for a first order phase transition that is stronger than in other 
treatments of the NJL model, is crucial to attain this result. 

We have studied the baryon number susceptibility and the specific heat which are related 
with event-by-event fluctuations of /i^ or T in heavy- ion collisions. In the two and three- 
flavor NJL models, for xb, we conclude that the obtained critical exponents around the 
CEP in both models are consistent with the mean field values e = e' = 2/3. For the specific 
heat we obtain nontrivial exponents 1/2 < a < 2/3 around the CEP, indicating a crossover 



of different universality classes 13|, |4^. This effect is more clearly visible for the critical 



exponent of the specific heat in the SU(2) version of the NJL model, where a crossover from 
a to «! is also observed. Nevertheless we notice that the values of a in the TCP and in the 
CEP are consistent within both versions of the NJL model. A better insight to the difficult 
task of the analysis of the phase diagram of QCD can be provided by an extension of the 
NJL model where quarks interact with the temporal gluon field represented by the Polyakov 



loop dynamics 
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5l| . Work in this direction is in progress. 



Concerning the behavior of the chiral partners in the vicinity of the CEP, we verified that 
the two Mott lines, respectively, for a and vr*^ cross at a point just bellow the CEP. On the 
other hand, there is a sharp decrease (increase) of the sigma (pion) meson masses at the 
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CEP, which reflects the nature of the second order phase transition at this point. 
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APPENDIX A 

In this appendix we give the integrals appearing in the meson propagators, in the vacuum 
and at flnite temperature and density, as well as some useful expressions concerning the a 
meson. 

The integrals I{ and II^Pq) are given by 

/i(r,/x.) = ^ J^{l-n.-n.), (Al) 

mPo, T, /X.) = / ^ W=l^' il-n.-n.), (A2) 

where Ei = + Mf is the quark energy. To regularize the integrals we introduce the 
3-dimensional cutoff parameter A. When Pq > 2Mj it is necessary to take into account the 
imaginary part of the second integral. It may be found, with help of the ie -prescription 
Pq Pq — ie. Using 

lim — — = V- + m5iy) (A3) 
t^o+y-ie y 

we obtain the integral 

Concerning the calculation of the propagator for the a meson, the projector Sab and 
the polarization operator Ilf^ matrices, in the case (quQu) = {QdQd), have the nonvanishing 
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elements 



5'33 = 9s- 9d {qs qs) , (A5) 
2 

'S'oo = 9s + ^9d {{qu qu) + {qd qd) + {qs qs)) , (A6) 

5'88 =9s- \^9d (2 {qu qu) + 2 {qd qd) - {qs qs)) , (A7) 

-^os = -^so = ^^^9d {{ququ) + {qdqd) - 2 (g^g^)) . (A8) 



Analogously, we get 



where 



n^o(^o) = I [Jim + JUPo) + 4(Po)] , (A9) 

nf8(Po) = \ [jL{p^) + J!d{Po) + 4Jf,(Po)] , (Alo) 

n^8(Po) = nfo(Po) = ^ [JL{Po) + ^|.(Po) - 2Jf,(Po)] , (All) 



4(Po) = 4[2Ji + [P2 - 4M^]/^^(Po)]. (A12) 
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